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MTIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

EXPERIMENTAL FLUTTER INVESTIGATION OF A THIN UNSWEPT 
WING AT TRANSONIC SPEEDS 


By George L. Pratt 


SUMMARY 


Zero-lift exj)erim.ental flutter data have been obtained in the Mach 
number range from 0.75 to 1.25 on a 1 . 5 -percent- thick, untapered, unswept, 
solid steel wing having hexagonal airfoil sections and wing len^h-to- 
chord ratios of 1.6l and 1.73. The wing was tested as a cantilever with 
and without a half -body-of -revolution fuselage. The test data indicated 
that reflected fuselage bow waves had an effect on the flutter results at 
Mach uTjmbers from approximately 1.05 to 1.20; hence, the discussion and 
conclusions are confined mainly to the data obtained with the fuselage 
removed. 

The results of the test of the wing at a length-to- chord ratio of 
1 . 6 l and with the fuselage removed indicated that calculated flutter 
speeds based on two-dimensional, incompressible-flow aerodynamic coef- 
ficients were approximately 15 to 20 percent less than the ejqjerimental 
values at Mach numbers from O .875 to 1.175; this difference increased at 
the higher test Mach numbers. An aspect-ratio correction resulted in cal- 
culated flutter speeds which were 10 to 15 percent greater than the experi- 
mental values at Mach n\mibers from O .875 to 1.175; this difference 
decreased at the higher test Mach nimbers. Except at the higher Mach num- 
bers, the analytic solutions satisfactorily predicted the variation of 
flutter speed with mass ratio, but the calculated flutter speeds based on 
two-dimensional aerodynamic coefficients were approximately I 5 to 20 per- 
cent less than the experimental values. 


INTRODUCTION 


Only a limited amount of research has been done experimentally and 
analytically to determine the effect of various design parameters on the 
flutter characteristics of wings at transonic speeds. (Some of the more 
recent experimental research is contained in refs. 1 to L.) The natirre 
of the airloads on wings at transonic speeds con5)licates flutter analysis 
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and increases the need for experimental flutter research. Such research 
provides the flutter data necessary to check various analytic methods of 
flutter prediction and provides data for eirpirlcal studies. An investi- 
gation has heen inaiigurated in the Langley 9- hy 12-inch hlowdown tunnel 
to provide experimental data indicating the effect of various wing geo- 
metric, elastic, and mass parameters on the flutter characteristics of 
wings at transonic speeds. 

The present paper contains the results of the part of this investi- 
gation concerning the zero-lift experimental flutter data obtained in the 
Mach number range from 0.75 to 1.25 on a 1.5-percent- thick, unswept, 
untapered, solid steel wing having hexagonal airfoil sections. Tests 
were made with the wing cantilever-mounted from the tunnel floor at a 
wing length-to- chord ratio of 1.6l. For some of the tests, the wing was 
cantilever-mounted from a half -body- of -re volution fuselage on the floor 
of the tunnel at aspect ratios of 3.75 and 4.00 which correspond to 
exposed-wing length-to- chord ratios of 1.6l and 1.75? respectively. 

The experimental flutter speeds, reduced flutter speeds, and flut- 
ter frequencies have been normalized by using flutter values calculated 
by the strip-analysis proced\are of reference 5 with aerodynamic coeffi- 
cients for incompressible flow. In addition, three-dimensional aerody- 
namic effects are taken into account for the case of the wing without the 
fuselage by applying an aspect-ratio correction determined by the method 
of reference 6. 

The present data sure the Initial flutter data obtained in the test 
facility, and because of the large ratio of wing plan-form area to test- 
section cross-sectional area, some comment on the reliability of the flut- 
ter data obtained is warranted. Althotigh some effect of the reflection 
of the model fuselage bow wave from the tunnel walls at certain Mach num- 
bers is Indicated in the present paper, it is believed that the test facil- 
ity is satisfactory for flutter-testing, particiilarly when the effects of 
the fuselage bow wave have been eliminated. This belief is based on the 
satisfactoiry results obtained in this facility in static tests of models 
of similar size and on the restilts of reference 1 which indicate that 
satisfactory flutter data can be obtained in transonic slotted wind 
tunnels . 


NOTATION 


A aspect ratio including body intercept 

a nondlmensional elastic-axis position measin-ed from mid- 

chord, positive rearward, semichords 
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b 


C 





wing semi chord, in. or ft (as required by dimensional 
considerations ) 

wing chord, 2b, in. 

first bending natural frequency, cps 
second bending natural frequency, cps 

uncoupled first torsion frequency about elastic eixis, cps 

structural dan5)ing coefficient 

structural daiirping coefficient for bending mode 


g 


a 


a 


^xx 

I 


M 

m 


structural damping coefficient for torsion mode 

wing mass moment of inertia about elastic axis per unit 
length, lb-sec^ 

wing section area moment of inertia, in.^ 

length of exposed wing measured perpendicular to wing 
root, in. 

Mach number 

mass of wing per unit length, Ib-sec^/sq in. 




dynamic pressure at flutter, ipV 

2 e ^ 


Ib/sq ft 


a 


t 

V 

X 

a 


nondimens ional radius of gyration of wing about elastic 


axis 




fraction of semichord 


wing thickness, in. 

stream velocity at flutter, ft/sec 

nondimens ional wing section center- of- gravity location 
measured from elastic axis, positive rearwaord of elas- 
tic axis, percent of semi chord 


CONFIDENTIAL 


4 CONFIDEOTIAL nacA RM L55AI8 

p air density at flutter, slugs/cu ft or Ib-sec^/in,^ (as 

required by dimensional considerations) 

p mass ratio at flutter, m/«pb2 

CD flutter frequency, radians/sec 

^ first bending natural frequency, radians/sec 

uncoupled first torsion frequency about elastic axis, 
radlans/sec 


Subscripts : 

e experimental value at flutter 

R based on strip-analysis method of reference 5 with 

two-dimensional aerodynamic coefficients 

A based on aspect-ratio correction of reference 6 


MODEL, APPARAITUS, AND TESTS 
Model 


The wing was mtapered, unswept, and of solid steel construction 
with a hexagonal airfoil section. The chord was 2 Inches and the 
thickness-to-chord ratio was O.OI5. Various other wing geometric and 
mass parameters are listed in table 1. Figure 1 shows a sketch of the 
model arrangement with the fuselage Installed. The fuselage was a hailf- 
body of revolution mounted on the floor of the tunnel with a l/l-inch 
shim to extend the model beyond the tvinnel- floor boundary layer. The 
wing was clamped rigidly at the wing- fuselage juncture and the aspect 
ratio could be changed by merely loosening the clamp and sliding the 
untapered wing through the fuselage. A sketch of the fuselage is pre- 
sented in figure 1. With the fuselage installed, the wing was tested at 
aspect ratios of 5.75 and l+.OO, These aspect ratios correspond to 
exposed-wing length- to-chord ratios z/c of I.6I and 1.73, respectively. 
The wing was tested without the fuselage by cantilever-moimting the wing 
from the tunnel floor at z/c=1.6l. 

The measured natural bending and torsion frequencies for the wings 
at the two aspect ratios and without the fuselage are presented in 
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tatle II, Inasmucli as the elastic axis and center of gravity coincide 
on the wing, these frequencies correspond to the xincoupled frequencies. 


Wind Tunnel 

The wind tixnnel used in this investigation was the Langley 9- hy 
12- inch blowdown txjnnel incorporating a transonic test section. The 
transonic test section is 7 inches high and 10 inches wide and is slotted 
longitudinally along the top and side walls. Without a model installed 
in the test section, tunnel sxirveys indicate that the maximum deviation 
of the Mach n\miber from the average test- section Mach momber over the 
length of the test section is i0.005 at M = 0.75 increases to ±0.020 
at M = 1 . 25 . Semi span models are cantilever-moxinted from the floor of 
the test section. 

The tunnel stagnation pressure can be varied from atmospheric pres- 
sure to approximately 51 Ib/sq in. abs by means of a throttling valve 
located upstream of the test section. The Mach number is controlled by 
means of a cylindrical plxinger located in the closed-wall part of the 
tunnel downstream of the test section. The plunger can be extended into 
the airstream and chokes the tunnel at the desired Mach number. This 
arrangement of throttling valve and Mach number control plunger permits 
the Mach number and stagnation pressure to be changed Independently. 

An air dryer and heater is Installed upstream of the test section 
to eliminate any condensation effects. 

In order to stop the flutter before the wing fluttered to destruc- 
tion, a quick-acting butterfly valve located upstream of the test section 
could be closed. The butterfly valve stopped the flow of air in the tun- 
nel in less than 1 second. 


Apparatus and Instrumentation 

Strain gages were installed on the surface of the wing near the 
root to indicate the bending and torsion frequencies and to establish 
the occttrrence of flutter. (Visual observation of the wing was restricted 
and the bending and torsion modes during flutter were difficult to 
discern.) The signal from the strain gages was acrplified and fed into a 
recording oscillograph to obtain the time histories of the flutter motion. 

The tunnel test conditions were obtained by photographically 
recording the ttmnel stagnation and static pressures (indicated on a 
manometer board) and the tunnel stagnation ten 5 )erature simultaneously 
with each oscillograph record. 
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Test Procediore 

The wing flutter points were obtained by approaching the flutter 
condition by two procedures: (l) setting the tunnel stagnation pressure 

at a desired value and changing the Mach number in small increments until 
the model fluttered and (2) setting the Mach number control plunger in a 
desired position and varying the test-section density by increasing the 
ttmnel stagnation pressure in small increments until flutter was obtained. 
In many instances it was found convenient to vary alternately both the 
Mach number and stagnation pressure dioring the course of obtaining a flut- 
ter point. 

During a test, the vibrations of the wing caused by the turbiilence 
of the flow in the tunnel were believed to be of sufficient amplitude so 
that external excitation of the wing to initiate flutter was not consid- 
ered necessary. 


RESULTS AND DISCUSSION 
Presentation of Data 


The flutter data obtained with the wing are presented in the fora 
of plots of several experimental flutter parameters (q., ^eA“e) 

plotted against the Mach number at which the data were obtained (fig. 2). 
for the purpose of providing a basis for comparing the results of the 
present investigation with results for other wings, an attempt has been 
made to account for the effects of such flutter parameters as 

M a, x_, and r^^ by normalizing the data by using calculated 
reference flutter speeds and frequencies. These normalized data are pre- 
sented in figure 3 as a function of Mach nimiber. In order to provide a 
direct comparison indicating the effect of the fuselage on the flutter 
characteristics of the wing, the values of Ve/Vp obtained for the wing 
with and without the fuselage at l/c = 1.6l are presented in figure k. 
Sections of the oscillograph-record bending and torsion traces from which 
the flutter frequencies were determined are shown in figure 5 ^"0^ several 
Mach numbers. In figures 6 and 7 the effect of the mass ratio on the 
wing flutter characteristics is indicated. (The quantity 

constant for any particular length-chord ratio and a, x^, and r^ were 

constant for the test wing throughout the present investigation.) 

Table III presents the experimental and analytical results obtained with 
the fuselage removed and the wing cantllever-moiinted from the tunnel 
floor at l/c = 1.61. 
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Calculated Flutter Speeds and Frequencies 

Two analytic methods have been utilized to calculate the reference 
flutter speeds and frequencies used to normalize the flutter data. The 
first method of analysis is based on two-dimensional, incompressible- 
flow, oscillating-alrfoll aerodynamic coefficients. The solution was 
obtained by applying the strip-analysis method of reference 5 tor a wing 
of zero sweep oscillating in two degrees of freedom and by employing 
the uncoupled first bending and uncoupled first torsion mode shapes of 
a uniform cantilever beam. The calculated reference flutter speeds and 
frequencies of this analysis eire denoted by the subscript R throughout 
the paper. 

The second analytic solution differs from the above solution only 
in that it employs the aspect-ratio corrections of reference 6. The 
reference flutter speeds and frequencies determined by this method of 
solution are denoted by the subscript A. 

In these two methods, the flutter determinant was solved at various 
values of reduced flutter speed V/bcu and mass ratio p. for the damping 
coefficient g. The faired vailue of V/bo) from a plot of V/bo) against 
g for a specific value of p at which g^ = g^ = g = 0 was taken to be 

the reference flutter point. 


Effect of Fuselage 

In the initial phase of the investigation, the wing was cantilever- 
mounted from a half -body of revolution as indicated in figure 1 . The 
flutter data obtained through the test Mach number range for aspect 
ratios of 4.00 (z/c = 1.75^ figs. 2(a) and 3(a)) and 5«75 (z/c = I.6I, 
figs. 2(b) and 5(b)) with this arrangement indicate a sharp deviation 
from the trend of the flutter points with Mach number in the proximity 
of M = 1.05 as shown by the large scatter in the data. A reason for 
this scatter in the data appeared to be that the reflection of the fuse- 
lage bow shock wave from the tunnel walls back to the model might be of 
siifficient strength to affect the flutter characteristics of the wing 
even at the lower supersonic Mach numbers. 

In order to eliminate any effects of the fuselage bow wave, the 
fuselage was removed and the wing was cantilever-mounted from the floor 
of the tunnel at the same exposed-wing length-to- chord ratio (z/c = I.61) 
as the A = 5*75 wing- fuselage configuration. The tunnel- floor bomdary 
layer at the root of the model was considered to have relatively little 
effect on the flutter characteristics of the model. 

With the fuselage removed, the deviation in the trend of the data 
with Mach number was eliminated near M = I.05. In addition, there was 
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a variation in the magnitude of the data from approximately M = 1.05 "to 
M = 1.20. Simple shock-reflection theory indicates that, in this Mach 
number range, shock waves from the fuselage bow might be expected to 
reflect back onto the wing from the tunnel walls. Figure k shows the 
TnaY^mnm variation in the flutter- speed ratio restating from fuselage 
effects to be approximately 0.08 at M = 1.17. Inasmuch as the fuselage- 
bow-wave reflections appear to have an effect on the flutter character- 
istics of the wing, the remainder of the discussion is confined to the 
data obtained with the fuselage removed and with the wing cantilever- 
mounted from the tunnel floor at an exp)Osed-wing length-to- chord ratio 
Z/c of 1.6l. 


Variation of Wing Flutter Characteristics With Mach Number 

An examination of the experimental flutter data for the wing at 
z/c = 1.6l (fig. 2(c)) indicates that less dynamic pressure q is 
required to flutter the model as sonic velocity is approached but that 
increasingly higher dynamic pressures are required before flutter will 
occur at Mach numbers greater than M = 1.175. Figure 2(c) also shows 
the variation of the ratio of the flutter frequency to the first natural 
torsion frequency uje/Hx variation of the reduced flutter 

speed VgA>^e number for the test wing and indicates a con- 

siderable effect of Mach number on the wing flutter frequencies. Fig- 
ure 5 presents sections of the oscillograph-record bending and torsion 
traces obtained at various Mach numbers from which the flutter frequen- 
cies were determined. 


Results of Analytic Solutions 

The variations with Mach number shown by the normalized flutter data 
(fig. 5(c)) indicate that there is a considerable effect of Mach number 
on the flutter characteristics of the wing. Some variation might be 
expected, particiaarly through the transonic speed range where the flow 
over the wing is in the transition range from subsonic to supersonic 
speed, because the attendant changes in the magnitudes and phases of the 
air forces and in the location of the aerodynamic center with these flow 
changes are not accounted for by the two-dimensional incongjresslble-flow 
aerodynamic coefficients . 

In general, the calculated flutter speeds based on two-dimensional, 
lncon 5 )ressible-flow aerodynamic coefficients Vr were approximately 
15 to 20 percent less than the experimental values from M = O .875 to 
M = 1.175 (fig. 3(c)). The aspect-ratio correction resiated in calcula- 
ted flutter speeds V/^ which were 10 to I 5 percent greater than the 
experimental values from M = O .875 ^ ~ 1»175 (fig. 3(c)). At the 
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Mach numbers above M = 1.175 "the flutter- speed ratios Vg/Vp and Vg/Vj^ 
increased with an increase in Mach number. 

The analytic solutions utilized did not accurately predict the flut- 
ter frequencies as indicated by the values of flutter-frequency ratio 
and (fig* 3(c)). 


Variation of Wing Flutter Characteristics With Mass Ratio 

In figures 6 and 7> various experimental and normalized flutter par*am- 
eters are plotted against the square root of the mass ratio at flutter. 
Since the effect of Mach number cannot be disregaxded in a cong)arison 
such as this, the data have been separated into groups corresponding to 
various parts of the Mach number range by using approximately the inflec- 
tion points of the Ve/VR plots of figure 5(c) to separate the various 
groupings (M < 0.975, 0.975 < M < I.050, I.050 < M < 1.175, and 
M > 1.175) . 

The data indicate that, in general, the variation of Vgydxi^ with 

was essentially the same as predicted hy the two analytic solutions 
except at the highest Mach numbers, but there was a difference in magni- 
tude (approximately I5 to 20 percent) between the experimental and cal- 
culated values based on the two-dimensional aerodynamic coefficients 
(fig. 6). The parameters V^/bo^ and cd^/o)^ differed both in magnitude 

and trend with from those predicted by the analytic solutions. 


CONCLUSIONS 


The results of the flutter tests of a thin unswept wing at transonic 
speeds indicated: 

1. For the wing cantilever-mounted without the fuselage at a length- 
to-chord ratio of I. 6 I, calculated flutter speeds based on incompressible- 
flow, two-dimensional aerodynamic coefficients were approximately I 5 to 

20 percent less than the experimental values at Mach numbers from 0.875 
to 1.175« This difference increased at the higher Mach numbers at which 
flutter data were obtained. 

2. An aspect-ratio correction resulted in calculated flutter speeds 
which were 10 to I 5 percent greater than the exp)erimental values at Mach 
n-umbers from 0.875 to 1.175. This difference decreased at the higher 
Mach numbers at which flutter data were obtained. 
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5. Except at the higher Mach numbers, the analytic solutions 
satisfactorily predicted the variation of flutter speed with mass ratio 
but the calciolated flutter speeds based on the two-dimensional aerody- 
namic coefficients were approximately I5 to 20 percent less than the 
experimental values. 

4. There was some effect of reflections of the fuselage bow wave 
from the timnel walls between Mach numbers of I.05 to 1.20. The maximum 
variation in flutter-speed ratio was O.O8 at a Mach nimiber of I.17O 
between the fuselage-on and fuselage-off conditions. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 15, 1955 • 
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TABLE I.- GEOMETRIC, MASS, AND INERTIA 
CHARACTERISTICS OF THE WING 


t. In 0.030 

c, in 1.998 

a 0 

0 

m, Ib-sec^/sq in 37 X 10“ 

or p)j,, Ib-sec^/in 11.9^ 

rtb , ^ 

I , in.^. 3.52 X 10“° 

xx' , 

I , lb- sec ^ 9*50 X 10"° 

a’ 

r 2 0.248 

a 


TABLE II.- WING NATURAL FREQUENCIES 


Z/c 

A 



fa 

(JJ 

a 

Fuselage 

1.61 

5.75 

89 

— 

552 

0.253 

On 

1.61 

— 

89 

— 

552 

.253 

Off 

1.73 

4,00 

74 

k6o 

326 

CVI 

CVJ 

• 

On 
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TABLE III.- EXPERItffiNTAL AND ANALYTICAL RESULTS FOR WING 
CANTILEVER-MOUNTED FROM TUNNEL FLOOR 


[z/c = 1.61] 


M 

IB/sq. 


Ve, 

ft/sec 

fe, 

qps 


Ve/t^e 


VeAE 


O.87J4. 

1,1*19 

79.8 

946 

153 

0.453 

11.80 

0.887 

1.205 

5.11 

.887 

1,411 

83.5 

964 

148 

.419 

12.45 

.862 

1.202 

5.22 

.907 

1,555 

89.9 

980 

150 

.425 

12.47 

.880 

1.178 

5.30 

.911 

1,381* 

90.7 

995 

145 

.406 

13.29 

.841 

1.194 

5.40 

.918 

1,31*6 

92.7 

992 

139 

.39^ 

15.64 

.818- 

1.177 

5.57 

.928 

1,525 

94.0 

991 

139 

.394 

15.62 

.819 

1.168 

5.57 

.9^0 

1,293 

100.7 

1,014 

156 

.385 

14.25 

.805 

1.154 

5.^8 

.9^4 

1,281 

100.2 

1,006 

151 

.371 

14.67 

.775 

1.149 

5.44 

.950 

1,276 

104.1 

1,024 

131 

.371 

14.94 

.778 

1.148 

5.54 

.960 

1,251* 

106.7 

1,028 

129 

.366 

15.22 

.769 

l.l4o 

5.57 

.966 

1,21*3 

114.0 

1,061 

126 

.358 

16.08 

.756 

l.l4o 

5.76 

.97‘^ 

1,255 

112.9 

1,049 

120 

.540 

16.69 

.718 

1.151 

5.68 

.983 

1,229 

115.9 

1,052 

ll4 

.325 

17.61 

.682 

1.128 

5.69 

.981^ 

1,257 

111.3 

1,051 

114 

.523 

17.61 

.681 

l.l4l 

5.69 

.988 

1,355 

107.^ 

1,064 

117 

.531 

17.56 

.696 

1.172 

5.75 

.990 

1,250 

113.7 

1,059 

115 

.526 

17.59 

.689 

1.139 

5.73 

.990 

1,263 

111.5 

1,054 

114 

.523 

17.67 

.681 

1.144 

5.71 

.995 

1,289 

111.0 

1,063 

116 

.329 

17. ^^9 

.694 

1.156 

5.75 

.997 

1,351 

106.8 

1,067 

115 

.529 

17.72 

.691 

1.192 

5.85 

1.000 

1,527 

107.9 

1,063 

119 

.337 

17.06 

.709 

1.171 

5.75 

1.011 

1,1*38 

101.4 

1,073 

122 

.546 

16.79 

.724 

1.219 

5.81 

1.012 

1,1*01 

106.1 

1,085 

120 

.541 

17.25 

.716 

1.207 

5.88 

1.012 

1,1*39 

101.4 

1,075 

121 


16.94 

.718 

1.220 

5.81 

1.012 

1,399 

104.8 

1,076 

120 

.540 

17.13 

.715 

1.202 

5.82 

1.022 

1,1*27 

105.1 

1,087 

122 

.546 

17.02 

.726 

1.215 

5.89 

1.029 

l,i*i*7 

104.4 

1,092 

121 

.545 

17.25 

.719 

1.222 

5.91 

1.055 

1,1*60 

104.8 

1,099 

125 

.549 

17.06 

.752 

1.229 

5.95 

1.056 

1,1*91 

102.1 

1,096 

122 

.546 

17.17 

.724 

1.242 

5.94 

1.056 

1,1*68 

104.4 

1,100 

122 

.5^^6 

17.21 

.726 

1.252 

5.96 

1.042 

1,1*57 

107.8 

1,115 

325 

.549 

17.28 

.734 

1.229 

6.05 

1.045 

1,478 

104.5 
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Figure 1.- 


Sketch of ving-fuselage 


arrangement. 
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(a) A = 4.00; z/c = 1.73; fliselage installed. 

Figure 2.- Experimental flutter characteristics plotted against Mach number. 
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(b) A = 5.75; z/c = I.6I; fuselage installed. 
Figure 2 .- Continued. 
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(c) z/c = 1.6l; no fuselage. 

Figure 2. - Concluded. 
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(a) A = 4.00; z/c = 1-75J fuselage installed. 

Figure 3-- Normalized flutter characteristics plotted against Mach nvimber. 
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(b) A = 5*75; z/c = 1.6l; fuselage installed. 

Figure 3* - Continued. 
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(c) l/c = 1.6lj no fuselage. 

Figure 5*- Concluded. 
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Figure 4.- Comparison of flutter speed ratios for wing with and without 
fuselage installed, z/c = 1.6lj A = 5.75. 
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Figure 5 . - Sections of oscillograph-record bending and torsion traces at 
several Mach nuirbers. z/c = I.6I; no fuselage. 
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Figure 6.- Effect of mass rafio on tlie fluffer characteristics. z/c = 1.6l^ 

no fuselage. 
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. - Effect of mass ratio on normalized flutter characteristics. 
Z/c = I.6I; no fuselage. 
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